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Abstract: An algorithm for generating a class of closed form 
solutions to the Navier-Stokes equations is suggested, with exam- 
ples. Of particular interest are those exact solutions that exhibit 
intermittency, tertiary Hopf bifurcations, flow reversal, and hys- 
teresis. 



1 INTRODUCTION 

The Navier-Stokes equations are notoriously difficult to solve. However, from 
the viewpoint of differential topology, the Navier-Stokes equations may be 
viewed as a statement of cohomology: the difference between two non-exact 
1-forms is exact. Abstractly, the idea is similar to the cohomology statement 
of the first law of thermodynamics. 



Q-W = dU (1) 

For the Navier-Stokes case, define the two inexact 1-forms in terms of the 
dissipative forces 



W D = f D » dr = p {iA7 2 V} • dr (2) 
and in terms of the advective forces of virtual work 



W v = fy • dr =p{dV/dt + grad(V • V/2) - V x curlV} • dr (3) 
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Then the abstract statement of cohomology, formulated as Wy — Wo = —dP , 
when divided by the common function, p, is precisely equivalent to an exterior 
differential system whose coefficients are the partial differential equations 
defined as the Navier-Stokes equations, 

{dV/dt + grad(V • V/2) - V x curlV} - {zA7 2 V} = -grad P/p (4) 

The cohomological constraint on the velocity field, V, is such that the kine- 
matically defined vector, f , 

f = fy - fo (5) 

is a vector field that satisfies the Frobenius integrability theorem [1]. That 
is, 

f • curl f = even though v« curl v 7^0. (6) 

The meaning of the Frobenius criteria is that the vector f has a representa- 
tion in terms of only two independent functions of {x,y,z,t}. The Navier- 
Stokes equations makes this statement obvious. One of these functions has 
a gradient, gradP, in the direction of the tangent vector to f, and the other 
function, p, is a renormalization, or better, a reparametrization factor for the 
dynamical system represented by f . 

These observations suggest that there must exist certain constraint rela- 
tionships on the functional forms that make up the components of any solu- 
tion vector field, V, (which usually does not satisfy the Frobenius condition 
in general) such that the covariant kinematic vector, f , is decomposable in 
terms of at most two functions. If such a constraint equation can be found in 
terms of the component functions that represent V, then its solutions may be 
easier to deduce than the direct solutions of the Navier-Stokes equations. For 
example, the constraint relation may involve only 1 partial differential equa- 
tion rather than 3. In fact such a single constraint relation can be found by 
imposing a type of symmetry condition on the system, a symmetry condition 
that expresses the existence of a two dimensional (functional) representation 
for the vector field, f . In this article attention will be focused on the two 
spatial variables, r and z, such that the solution examples will have a certain 
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degree of cylindrical symmetry. As these solutions involve dissipative terms 
with a kinematic viscosity coefficient, u, they are not necessarily equilibrium 
solutions of an isolated thermodynamic system. 

Closed form solutions are few in number [3], but it appears that many of 
the known steady-state solutions to the Navier-Stokes equations fall into the 
following class of systems: Consider a variety {x,y,z,t} with r 2 = x 2 + y 2 . 
Consider three arbitrary functions, Q(r,z) and $(r, z, i), and A(r,z) which 
are defined in terms of two independent variables spatial variables, (r, z), and 
time. Define the flow field, V in cylindrical coordinates as, 

V = A(r, z)u z + 0(r, z)u r + $(r, z, *)u^/r, (7) 

where is a unit vector in the azimuthal direction. Note that this vec- 
tor field does not necessarily satisfy the Frobenius theorem. Note that for 
simplicity, the only time dependence permitted is in the azimuthal direction. 

Substitution of this format for V into the equation for f will yield a vector 
equation of the form 

f = a(r, z)\i z + P(r, z)u r + 7(7-, z, t)u^/r. (8) 

The Pfaffian form W = f o dr will become an expression in two variables 
if the azimuthal factor 7(7", z, t) is constrained to the value zero. In other 
words, a single constraint on the functions, 9(r, z) and $(r, z, t), and A(r, z), 
defined by the equation 7(7-, z, t) = 0, can be used to reduce the Pfaffian form 
to the expression 

a(r, z)dz + (3(r, z)dr = -dP/p(r, z) (9) 

The left hand side represents a Pfaffian form in two variables, and therefore 
always admits an integrating factor. It is this idea that is used to find new 
solutions to the Navier-Stokes equations. First a solution to constraint equa- 
tion is determined. Then the Cartan 1-form of total work is computed. The 
1-form is either exact, or can be made exact by an appropriate integrating 
factor. If the 1-form is exact then the Pressure is obtained by integration. 
It the 1-form is not exact a suitable integrating factor is found, and that 
integrating factor represents a variable fluid density, p. For a given choice of 
integrating factor, the Pressure is again obtained by integration. 
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It is also useful to consider a rotating frame of reference defined by the 
equation 

n = uju z . (10) 

It is the choice of rotational axis that defines the cylindrical symmetry. For 
such rotating systems the same technique will insure that the flow field, V, 
is a solution of the Navier-Stokes equations in a rotating frame of reference, 

dV/dt + grad(Vo V/2) - V x curlV (11) 
= -gradP/p + z/V 2 V - 2Q x V - Q x (Q x r) (12) 

By direct substitution, into the Navier-Stokes equation above, of the pre- 
sumed format for the velocity field V yields an expression for 7(7", z) in terms 
of the three functions 0(r, z) and $(r, z), and A(r, z) : 

7 ( r ,z) = {d<5>/dt + A(r,z)d<$>/dz + Q(r,z)(d<$>/dr-2ujr) (13) 
-v{d 2 $/dz 2 + <9 2 $/<9r 2 - (d$/dr)/r}. (14) 

Similar evaluations of the standard formulas of vector calculus in terms of the 
assumed functional forms for the velocity field lead to the useful expressions: 

divV = d&/dr + 0/r + dA/dz (15) 

curl V = {d<S>/dr u z - d$/dz u r }/r + {80/ dz - dA/dr} u 

curlcurlV = {-d 2 A/dr 2 + d 2 Q/dzdr}u z + (16) 
{-d 2 Q/dz 2 + d 2 A/dzdr} u r + 
{-d 2 $>/dz 2 - d 2 $/dr 2 + (d^/dr/r)} u^/r 

VxcuHV = &{d&/dz-dA/dr}u z + (17) 
A{dA/dr- dQ/dz} u r + 

{(1/r 2 ) grad ($ 2 /2) - {Ad$/dz + 09$/9r} u^/r 
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grad(ym\)/2 = {QdQ/dz + AdA/dz + <^(d^/dz)/r 2 }u z + (18) 
{QdQ/dr + AdA/dr + $(<9$/<9r)/r 2 - $ 2 /r 3 }u r 



grad(divV) 



{d 2 @/dzdr + (dG/dz)/r + d 2 A/dz 2 }u z + 
{d0 2 /dr 2 + d 2 A/drdz + {dQ/dr)/r - Q/r 2 }u r 



(19) 



It is remarkable that many solutions to the Navier-Stokes equations then 
can be found by using the following algorithm: Choose a functional form 
$(r,z) of interest and then deduce functions A(r,z) and Q(r,z) to satisfy 
the azimuthal constraint, 



The flow field V so obtained is therefore a candidate solution to the 
compressible, viscous, three dimensional Navier-Stokes equations for a system 
with a density distribution, p and a pressure, P. The components of flow field 
so determined then permit the evaluation of the coefficients of the Pfafnan 
form 



If the expression is not a perfect differential, then use the standard methods 
of ordinary differential equations to find an integrating factor, p(r,z). The 
integrating factor represents the density distribution of the resulting Navier- 
Stokes solution. The Pressure follows by integration. 

This method is demonstrated in the next section for the known viscous 
vortex examples reported in Lugt. In addition, several new closed form exact 
solutions are generated by the technique. Among these closed form solutions 
are exact solutions to the Navier Stokes equations (in a rotating frame of 
reference) that exhibit the bifurcation classifications for N = 3 as given by 
Langford [2]. In particular, exact, non-truncated solutions are given that 
represent the trans-critical Hopf bifurcation, the saddle-node Hopf bifurca- 
tion, and the hysteresis Hopf bifurcation. It has been long suspect that many 



l(r,z,t) = 0. 



(20) 



W = a(r, z)dz + /3(r, z)dr 



(21) 
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phenomena in hydrodynamics exhibit Hopf bifurcation; now these exact so- 
lutions to the Navier-Stokes equation formally justify this position, and are 
especially interesting for the understanding of slightly perturbed Poiseuille 
flow and the onset of turbulence in a pipe. 

2 EXAMPLES 

In the following examples, the vector field specified has been used to compute 
the various terms in the Navier-Stokes equations. The algebra has been 
simplified by use of a symbolic computation program written in the Maple 
syntax. For each example, the two vector components that make up the work 
one form have been evaluated and are displayed with the solution. For the 
divergence free cases, the pressure function also has been computed. First, 
known solutions are exhibited, and are shown to be derived from the above 
technique. Then a few new solutions are exhibited. 

2.1 Old solutions 

2.1.1 Example 1. The Rankine Vortex 

$(r,z) = a+ (b + tu)r 2 , 9 = 0, A = l (22) 



f v = {-( a + 6r 2 ) 2 /r 3 }u r + {O}u /r + {0}u, (23) 



f D = {0}u r + {O}u /r + {OK (24) 

This flow is a solution independent of the kinematic viscosity coefficient 
(the velocity field is harmonic, as = 0) and therefore could be construed 
as an equilibrium solution. This solution, for a and b equal to piecewise 
constants, will generate the Rankine vortex. 

As the flow is isochoric (div V = 0), the steady pressure can be determined 
by quadrature, and is given by the expression, 

P = 1/2(&V + Aabr 2 ln(r) - a 2 )/r 2 (25) 
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2.1.2 Example 2. Diffusion Balancing Advection. 



$(r, z) 



a + br 2+m/u , 0(r, z) = m/r, A = 1, w 







(26) 




(a + b^ 2u+m ^ v ) 2 /r 3 }u r + 
]/u m{2v + m)/z/r 2 }u^/r + {0}u 2 



(27) 



f D = {0}u r + {br i2u+m)/u m{2u + m)/ur 2 }u^/r + {0}u z 



(28) 



In this case the Laplacian of the vector field is not zero, but the dissipative 
parts exactly cancel the advective parts in the coefficient of the azimuthal 
field, thereby satisfying the constraint condition. As the functions depend 
only on r, the integrability (gradient) condition is satisfied, and these solu- 
tions obey the Navier-Stokes equations for a system of constant density. The 
Pressure function may be computed as 



The solutions are cataloged in Lugt. As these solutions explicitly involve 
the kinematic viscosity, u, they cannot be equilibrium solutions to isolated 
systems. Instead they represent steady state solutions, far from equilibrium. 
A special case exists for m/u = —2. 

2.1.3 Example 3. Burger's Solution, but with Helicity and Zero 
Divergence. 



$(r,z) = fc(l -e- ar2/2u ), G(r,z) = -ar, A = U + 2az, u = (30) 

f v = {-{ke { - ar2/2v) +r 2 a-k)(ke { ' ar2 / 2 ^ - r 2 a - k)/r 3 }u r + 
{kra 2 /u e { ~ 1/2ar2/u) }u^/r + 
{2{U + 2az)a}u z 



P = 



(-ub(4a(m + + bmr (2u+m)/u) )r (2u+m)/u) 
(m(u + m)(a 2 + m 2 ))/{2m(v + m)r 2 ) 



(29) 
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f D = {0}u r + {-kra 2 /v e^ 2ar2 ^}u^/r + {0}u z (31) 

This solution corresponds to a modification of Burger's solution and ex- 
hibits a 3- dimensional flow (in 2-variables) in which the diffusion is balanced 
by convection to give azimuthal cancellation. The Burgers solutions has 
been modified to exhibit zero divergence. This flow in a non-rotating frame 
of reference exhibits a helicity. 

Helicity = (U + 2az) (ka/u)e^ 1/2ar ' '' /v) (32) 
2.2 New Solutions 

2.2.1 Example 4. A Beltrami Type Solution 

$(r, z) = r 2 cos(z/a), 0(r, z) — r sm(z/a), A(r, z) = 2a cos(z/a), oo = 

(33) 

fy = {— r}u r + {Oju^/r + {—Aacos(z/a)sin(z/a)}u z (34) 

fo — v /a 2 \{—r sin(z / a)}u r + {—rcos{z/a)}vi ( f ) /r + {— 2a cos(z/a)}u z ] 

(35) 

This solution is a Beltrami-like solution, has zero divergence, and can 
be made time harmonic by multiplying the velocity field by any function of 
t. The flow exhibits Eckman pumping and has a superficial resemblance 
to a hurricane. The time independent steady flow is a strictly Beltrami 
(curl v = a v)with the vorticity proportional to the velocity field. In all 
cases the helicity is given by the expression, 

Helicity := (r 2 + 4a 2 cos(z/a) 2 ) /a. (36) 

The kinetic energy is a/2 times the helicity, which is a times the enstrophy. 
The Pressure generated from the Navier Stokes equation is given by the 
expression 

P = l/2(r 2 + [r 2 {yja 2 ) - Av)sin(z/a) + 4a 2 sin(^/a) 2 ) (37) 
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2.2.2 Example 5. A Saddle Node Hopf Solution 



${r,z)=ur 2 , Q(r,z) = r(a + bz), A(r, z) = U - dr 2 + Bz 2 (38) 

The components of the advective force and dissipative force are given by 
the expressions, 

f v = { r (a + bz) 2 + (U - dr 2 + Bz 2 )rb}u r + (39) 
{0}u^/r + 

{-2r 2 (a + bz)d + 2(U - dr 2 + Bz 2 )Bz}u z 

and 

f D = {0}u r + {0}u^/r + u{-4d + 2B]}u z (40) 
The divergence of the velocity field is given by the expression: 

divV :=2{a + (b + B)z} (41) 
The helicity of the flow depends upon the rotation, u, 

Helicity : u(+r 2 b + 2U- 2bz 2 ) (42) 

but remarkably changes for finite values of r and z, depending on mean flow 
speed, U. 

Note that when b = 0, B = 0, a = 0, the solution is equivalent to the 
standard incompressible Poiseuille solution for flow down a pipe. The vector 
velocity field is not harmonic, but vector Laplacian of the velocity field is a 
constant. 

Without these constraints, it is remarkable that the ordinary differential 
equations that represent the components of the velocity field are in one to one 
correspondence with the saddle node - Hopf bifurcation of Langford. That 
is, the ODE,s representing the Langford format for the SN-Hopf are given 
by the expressions: 
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dz/dt = A(r, z) = U - dr 2 + Bz 2 (43) 
dr/dt = 0(r, z) = r(a + bz) 
d9/dt = uj 



.This first order system which exhibits tertiary bifurcation is associated 
with an exact solution of the Navier Stokes partial differential system in a ro- 
tating frame of reference. In principle, the method also relaxes the constraint 
on incompressibility, and allows a density distribution, or integrating factor, 
to be computed for an exact solution to the compressible Navier-Stokes equa- 
tions which can be put into correspondence with saddle node-Hopf bifurcation 
process. 

This example exhibits isochoric (divV = 0) flow for B + b = 0, a = 0. The 
steady isochoric pressure is then determined by quadrature, and is given by 
the expression, 



where the constant U can be interpreted as the mean flow down the pipe. Part 
of the pressure is due to geometry, and part is due to the kinematic viscosity. 
Note that the pressure is independent from the viscosity coefficient when the 
velocity field is harmonic; e.g. when (2d — b) = 0. As the vector Laplacian of 
the velocity field determines the dissipation in the system, intuition would 
say that the harmonic solution is some form of a limit set for the otherwise 
viscous flow. 

2.2.3 Example 6. A Transcritical Hopf Bifurcation 



$(r, z)=cor 2 , Q(r,z)=r(A-a + cz), A(r, z) = br 2 + Az + Bz 2 (45) 



P = b(dr 4 /2 - (r 2 - 2z 2 )U - bz 4 )/2 - v(Ad - 2b)z 



(44) 



f v = { r (A - a + czf + (br 2 + Az + Bz 2 )rc}u r + 
{OK/r + 

{2r 2 (A - a + cz)b + (br 2 + Az + Bz 2 )(A + 2Bz)}u z 



(46) 
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fo = {0}u r + {0}u^/r + u{Ab + 2B}u z 



(47) 



This example exhibits isochoric (divV = 0) flow for a = 3A/2 and B = —c. 
The steady isochoric pressure is then determined by quadrature, and is given 
by the expression, 

P = -l/4c6r 4 - 1/84 V + l/2A 2 z 2 - Az 3 c + l/2cV - v{Ab - 2c) z. (48) 

Again it is apparent that the pressure splits into a viscous and a non- viscous 
component, and when the flow is harmonic (26 — c = 0), the pressure is 
independent from viscosity, and there is no dissipation in the flow. 

The transcritical Hopf bifurcation is represented by the Langford system 

dz/dt = A(r, z) = br 2 + Az + Bz 2 (49) 
dr/dt = 0(r, z) = r(A - a + cz) 
dO/dt = uj 

2.2 A Example 7. A Hysteritic Hopf Bifurcation 

$(r, z) = cur 2 , Q(r,z)=r(a + bz), A(r, z) = U - dr 2 + Az + Bz 3 (50) 



f v = {r(a + bzf + (U - dr 2 + Az + Az 3 )rb}u r + (51) 
{0}u^/r + 

{-2r 2 (a + 6^)rf + (U - dr 2 + Az + Az 3 )(A + 3Az 2 )}u z 



f D = {0}u r + {O}u /r + z/{-4rf + 6Az}u z (52) 

This system has the remarkable property that the vector Laplacian changes 
sign at a position z = 2d/3A down stream. There is no global way of making 
this solution isochoric, for the divergence is equal to 
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divV = (A + 2a) + 2bz + 3Az 2 . 



(53) 



The hysteretic Hopf bifurcation exhibits what has been called intermit- 
tency. The Langford system is 
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dz/dt 
dr/dt 
dO/dt 



A(r,z) 
Q(r,z) 



U-dr 2 + Az + Bz 3 
r(a + bz) 



(54) 
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